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We show that the spin dynamics of the 4-channel Kondo 
model is equivalent to the one of a model with a spin 1/2 
impurity coupled to spin 1 conduction electrons. By abelian 
bosonization of the latter model we find a kind of Toulouse 
limit analogous to that recently discussed by Emery and 
Kivelson for the 2-channel case. We analyse the model by 
exploiting the close analogy to the problem of an impurity in 
the Luttinger liquid. 
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The overscreened multi-channel Kondo model has 
been the subject of an intensive analysis after Nozieres 
and Blandin's discovery of its non-Fermi-liquid behav- 
ior pi. Various methods have been used to tackle this 
problem: Bethe ansatz technique O], large n (the num- 
ber of channels) expansion |jh,|| and conformal field the- 
ory P,pi. The main feature of the overscreened Kondo 
model is a non-analytic temperature and magnetic field 
dependence of the impurity correction to the free energy, 
which manifests in the fractional exponents characteriz- 
ing the power law behavior of, e.g., susceptibility and 
specific heat. The exact exponents and the asymptotic 
behavior of the dynamic correlation functions have been 
determined by means of the Bethe ansatz and conformal 
field theory. 

Recently, Emery and Kivelson found an alternative ap- 
proach to the problem based on simple abelian bosoniza- 
tion g . It is analogous to the well known Toulouse limit 
for the standard single channel Kondo model and 
allows one to re-obtain all the exact results plus addi- 
tional information on dynamic properties. However their 
method is only applicable to the case n = 2. 

In this paper we show that a similar approach can be 
devised for the four channel model. The corresponding 
Toulouse limit is not exactly solvable on the contrary to 
the n = 1,2 case. Nevertheless, we are able to determine 
its low energy properties, which are sufficient to iden- 
tify all the exponents which control the low temperature 
behavior of the thermodynamic quantities. 

The hamiltonian for the n-channel Kondo model, H = 
Ho + Hexc, consists of the term describing free electrons. 



and the exchange term, 
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where la are the exchange constants, Sa are the impu- 
rity spin-1/2 operators, and J°(x) are the electron spin 
currents (densities): 



r(x)=^4,(x)r,V^..'(^), 
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r"^ being spin 1/2 matrices. As usually (see, e.g., p), we 
have retained the s-wave scattering only, linearized the 
fermion spectrum, and replaced the outgoing and ingo- 
ing waves with n-copies of right-moving electron fields 
ipiaix) defined for — oo < a; < -|-cx). Following |q] we 
introduce charge and "fiavor" currents, 

J{x) =J2^l{x)^^^{x), J^{x) ^J2^l{x)tf^^j^,{x) 
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tf; being generators of SUin) group, and re-write the 
free part of the hamiltonian as a sum of three commuting 
terms: 
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This allows one to formulate the problem entirely in 
terms of the electron spin currents J{x). The informa- 
tion about the number of channels is contained in the 
commutation relations obeyed by these currents [which 
follow from their definition Eq.(n])]: 



[J^ix), rm = ie'""'S(x)J''(0) - i—5'"'6'(x) 
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indicating that J°-{x) form SU {2)n Kac-Moody algebra 
(see, e.g., Q). In order to study magnetic properties, we 
include also a magnetic field h, which gives rise to an 
additional term in the hamiltonian 
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Our idea is to model the commutation relations (0), 
for particular values of n, by some free fermion repre- 
sentation of the currents, which will eventually be more 
convenient to work with (in a sense to be specified later) 
then the original representation given by Eq. (nl) . 

The following currents: 



J+{x) = V2 [4{x)^_,{x) + ^Jli^)^^) 
r{x) = i^\{x)i;^(x) - ^li(x)^„i(a:). 
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where ^^(x), m — —1,0,1, are spinless right-moving free 
electron fields, can be shown to obey the commutation 
relations (g) with n — 4. One can think of Eq.(^ as a 
representation of spin 1 electron currents. 

Since the spin dynamics is completely determined by 
the commutation relations obeyed by the electron spin 
currents, we can equally well work with one or the other 
set of J"'s, Eq. (|l|) or (||); the spin properties should 
be the same. This has also an interesting outcome that 
the spin behavior of the four channel Kondo model is 
equivalent to the one of a spin 1/2 impurity coupled to 
spin 1 electrons. The latter problem is interesting by 
itself, being the simplest realization of a Kondo model 
with integer spin conduction electrons. 

What we are going to show is that the representation 
(^ has the advantage of having a simple Toulouse limit, 
which the original representation (|lj) seems not to ex- 
hibit. We bosonize the fermion fields according to 

where (t>m{x) are free right- moving Bose fields and a 
stands for a high-energy cutoff (lattice spacing). Com- 
bining these Bose fields as 

<j>{x) = [M^] + M^) + cl)~i{x)]/V3, 

0,(a;) = [0i(x)-0_i(x)]/y2, 

^f{x) = [01 (x) - 2Mx) + 0-i(a;)]/\/6, 

we see that the field (t>{x) decouples from the impurity, 
and the part of the hamiltonian, Hk, which is relevant 
for the Kondo effect, takes the form 
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where Ho{ip) is the free Bose field hamiltonian. 
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and Ix = ly = I± is assumed (for convenience everywhere 
we do not indicate explicitly the x-dependence of the 
Bose fields, we intend them at the impurity site a; = 0). 
The next step is to employ the canonical transforma- 
tion: 

Hk^UHkU\ C/ = exp{-z5,0,/\/2}. (5) 

The transformed hamiltonian reads: 

Hk = Ho{(f>s) + Ho{c^f) + ^^S^ cos | ^^(bf J + 



SH 



h 



N/27r 



dxdx4)s{.x) 



(6) 



5H - XSAc, 



(7) 



where v^ttA — Iz ~ t^vf ■ In order to get rid of the bulk 
contribution of the magnetic field, we shift the Bose field 
according to 

dx<i>s{x) -^ dx4is{x) - V2h/vF ■ 
This shift changes 
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and generates a constant ~l/2xoh? times the length of 
the system, where xo — I/t^vf is the uniform spin sus- 
ceptibility of free conduction electrons. 

We see that the line A = is indeed analogous to 
the Toulouse limit for the single-channel Kondo prob- 
lem and to the Emery-Kivelson line for the two-channel 
one. Along this line, the canonical transformation (|^) 
leads to the decoupling of the impurity degrees of free- 
dom from the conduction electron ones [the impurity spin 
component Sx commutes with the hamiltonian (o) and 
hence loses its dynamics]. However, the hamiltonian for 
the phase field (j)f remains nontrivial. Let us write it in 
the form: 



(9) 



H{(bf) = Ho{(bf) ± Ii_ cos (^v/2.90/ j /V2T:a 



where for convention the ± sign refers to the conserved 
spin component 5*2; equal ±1/2. In our case g — 3/4. 
For the two-channel problem, the same hamiltonian was 
found with 5 = 1/2 in Ref. [|. 

We observe that the hamiltonian (0) is equivalent to 
the one describing a backscattering of electrons from a 
single impurity in the spinless Luttinger Liquid. The 
hamiltonian for this model, expressed in terms of the 
original Fermi operators, is: 



H = Hi 
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where the field 'iB.(L) refers to right (left) moving 
fermions, and H^uik is the interacting Luttinger liquid 
hamiltonian. If one bosonizes, according to the standard 
rules, these Fermi fields and gets rid of the interaction by 
performing a Bogoliubov rotation, the resulting hamilto- 
nian is indeed Eq. (0) pcf| . The parameter g in that equa- 
tion is a measure of the strength of the fermion-fermion 
interaction, being <? = 1 for non-interacting fermions and 
g < l{g > 1) for repulsivc(attractive) interaction. The 
problem of an impurity in the Luttinger liquid recently 
became quite popular in the context of the physics of 
one-dimensional quantum wires, and it was studied in 
great detail [10| . We will make use of some known results 
for this problem to analyse the, as we showed, equivalent 
Kondo model. Since in the specific case we are consider- 
ing the parameter g = 3/4 < 1, the impurity scattering 



operator is relevant and I± flows to infinity under scal- 
ing transformation. The resulting (stable) fixed point 
describes a perfectly reflecting barrier, which cuts the 
system into two semi-infinite lines a and b. Exactly at 
the fixed point, that is for everything regarding zero en- 
ergy and temperature properties, the two regions a and 
b are disconnected. At any finite frequency there is a 
residual tunneling across the barrier, described by the 
operator 



O = *1(0)*,(0) + i/.c. 
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where ^a(&) is the electron field referred to region a{b). 
This operator scales to zero approaching the fixed point 
as cji/9-i, with correlation function (O(i)O(O)) -- i^^/^ 
(see Ref. ^). 

Coming back to the Kondo model (||) and (Q), we no- 
tice that on the line A = the total uniform susceptibility 
X coincides with xoj that one of the conduction electrons 
in the absence of the magnetic impurity. Thus the im- 
purity susceptibility defined as Ximp = X ~ Xo vanishes 
on this line (the same occurs at the equivalent line of 
two-channel Kondo model ||5|,|6[|). Similarly the impurity 
contribution to the total specific heat at A = vanishes. 
Therefore, in order to study the physical properties of 
the Kondo model we need to perform a perturbation ex- 
pansion in the correction 5H, Eq. (@), to the hamiltonian 
(^ around the solvable line A = Q. As a first step, 
we calculate the dimension x of the impurity spin 5'^, 
which determines the large time decay of the correlation 
function 



{T{S,{t)SM))^Ct- 
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where C is a numerical coefficient. 

Since Sz is a sum of lowering and raising operators for 
the states labeled by the Sx component of the spin, its ef- 
fect is simply to reverse the sign of the backscattering po- 
tential Ij^ in the hamiltonian (g) . Therefore the problem 
of finding the dimension of Sz is similar, but not exactly 
equivalent, to the problem of calculating the X-ray edge 
exponents in the Luttingcr liquid. The latter problem 
corresponds to a sudden switching of the backscattering 
potential on(ofF) and has recently been studied ||ll|]. The 
difference in our case is that, instead of switching on(off) 
the scattering potential, we must cope with the sudden 
change of its sign. For non-interacting electrons {g — 1) 
the dimension of this operator can easily be calculated 
by a standard approach based on phase-shift arguments 
| jl2[ , and it turns out to be a; = 1/4. For the value g = 1/2 
the exponent a; = 1/2 was recently found in Ref. [p|. For 
arbitrary g < 1, we find x = l/4(;. This result, although 
related to the x-ray edge singularity, differs in a striking 
way, specifically the maximum value of the exponent is 
not bounded by the scattering phase shift value in the 
unitary limit. 



The exponent x = 1/4(7 is calculated by noticing that 
the unitary operator W = exp{z7rJ4/2} changes sign of 
the backscattering potential provided that the operator 
Jt is the total fermionic current |l3] , which is defined as 
the difference between the total numbers of right and left 
moving fermions and satisfies 



[Jt,cl)f{x)]^i./2/^. 



Therefore 



{T {Szit)Szm) ^ (T iW{t)W{0))). 

Due to the backscattering potential, the current Jt is not 
a conserved quantity but acquires a dynamics which has 
been analysed in Ref. |14|. From this analysis, Jt turns 
out to be essentially a gaussian variable with correlation 
function behaving at large time as 

{T{Jtit)JtiO))) = i2/n^g)\nt. 

Hence the dimension of W (i.e. Sz) is indeed x = l/ig. 
The dimension of the operator SH, Eq. (M), is 1 -I- 1/4(7 > 
1, thus being irrelevant at the stable zero temperature 
fixed point X — 0, I± = oo. The finite temperature 
behavior of the system is described by the irrelevant op- 
erators which move the hamiltonian away from this fixed 
point, specifically the operator SH, which brings the sys- 
tem away from A == 0, and the operator XO [see Eq. (|lO[)], 
corresponding to a deviation from Ij_ = oo, where A is 
a non-universal parameter. As we showed, the former 
operator has dimension 1 -I- l/4f/ and the latter one has 
dimension 1/g. For our case, g = 3/4, the two dimen- 
sions coincide, which implies that both operators have to 
be considered on equal footing. The two-channel Kondo 
model corresponds to 5 = 1/2, then O is more irrelevant 
than SH and can be neglected. 

Armed with this results, we can evaluate the first non- 
vanishing correction to the free energy due to a finite A 
and A in the presence of a magnetic field. The contribu- 
tion proportional to A^ is 

6F^Z--\l dT{T[dH{T)SH{Q)])^-^J^ dr 
{TlSzir) IdxMr) - ^1 ^.(0) ( d^MO) - ^ )]>, 
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[3 being the inverse temperature T. Since <j)s and Sz are 
independent, the above correlation function factorizes, 
and, by substituting the imaginary time representation 
of the correlation function (|ll|) , we obtain 
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The term proportional to /i^, which defines the mag- 
netic susceptibihty, is convergent (a; = 1/3). The other 
term, whose temperature dependence relates to the spe- 
cific heat, is instead singular, and needs to be regularized 
by introducing an ultraviolet cut-ofi'. The cut-off de- 
pendent part is analytic in temperature (actually it is a 
constant apart from 0{T^) terms), so that the regulariza- 
tion procedure does not affect the leading non-analytic 
T dependence of the free energy. The result is 
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where the coefficient 
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Analogously the contribution to the free energy propor- 
tional to A^ is: 
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By substituting the value a; = 1/3 into the above expres- 
sions, we get the impurity contribution to the specific 
heat and susceptibility: 



L^« 



2^2/3 



/(.imp 



T-1/3, 



in agreement with exact critical exponents found by dif- 
ferent methods |,|,|. 

Another quantity of physical interest is the Wilson ra- 
tio, which can also be evaluated, and for which wc find 
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where the bulk value is taken as Txo/Co — 3/(477^). It 
is worthy to note that, as follows from our solution, the 
Wilson ratio is not universal (in spite of the spin degrees 
of freedom have been kept decoupled from the charge 
and flavor ones). Since A and A describe how one moves 
away from the fixed point in the Iz and I± directions 
respectively, the non universality of Rw simply reflects 
the spin anisotropy. In the spin isotropic case the value 
Rw = 8 was found in Ref . Isj . The relevance of the spin 
anisotropy in the four channel model is a consequence 
of the operators SH and O having the same dimension. 
This contrasts to what happens in the two channel model, 
where only the former operator contributes, and in the 
conventional single channel Kondo model [[l5| ; in both 
cases spin anisotropy does not affect the Wilson ratio. 

To conclude, we have shown that the Toulouse limit 
can be realized in the four channel Kondo model or, 
equivalently, in the model where spin 1 conduction elec- 
trons couple to a spin 1/2 impurity. 
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